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Oh! Abstract 



We study the mixing of elementary and composite particles. In quantum field theory the 
mixing of composite particles originates in the couplings of the constituent quarks and for 
^ ' neutrinos in self-energy diagrams. In the event that the incoming and outgoing neutrinos 

■ have different masses, the self-energy diagrams vanish because energy is not conserved but 

the finite decaying widths make the mixing possible. We can consider the neutrinos to 
be "fuzzy" states on their mass shell and the mixing is understood as the overlap of two 
wavef unctions. These considerations restrict the mass difference to be approximately equal 
to or smaller than the largest of the two widths: |Mj — Mj\ ^ max\Ti^ Fj}. 



1 Lagrangian with Majorana neutrinos. 



During the last few years a lot of attention was paid to the possibility of creating a baryon 
asymmetry through leptogenesis. The proposed schemes introduce heavy Majorana neutrinos 
with CP-violating couplings, where both the so-called "direct" and "indirect" contributions to 
leptogenesis were considered. In calculating the lepton asymmetry, the authors consider the 
Standard Model (SM) with the usual particle content plus three Majorana neutrinos, which are 
singlets under the weak SU(2)-group - . 



The part of the Lagrangian with Majorana neutrinos consists of the Majorana mass term and 
the Yukawa interactions of these neutrinos with leptons and Higgs bosons: 

^ = ^ E MMN, + KJ^a<pPRN^ + KMPl lLa(l>^ + h.C. (1) 
i a,i a,i 

In this Lagrangian are left-handed lepton doublets of the SM, cj) is the Higgs doublet of the 
SM, a, 2 = 1,2,3 denote the index of fermion generation and Ni is the self-conjugate Majorana 
field: 

One should remember here that {NmY is a left-handed antiparticle (iVa)"^ = {Ni)^. 
Another (more explicit) way of writing this Lagrangian is used in 0: 



+ E ^ai(iVffi)^(/La)'=0 + E K^ilLaYiNni, 



(3) 



One can easily show that the two Lagrangians (P and (|^) are identical. In fact, by the definition 
of charge-conjugate fields one can show that 



{lLaY<P^{NR{)' = -{iLaV C''<P+CNr,' = -{iLaV (f>^ Nr,' = Nr^+Il^. (4) 

Notice that the mass term in ([^) violates lepton number by two units and the Yukawa inter- 
action terms violate the CP-symmetry. In general, it is desirable to have a gauge theory with a 
symmetry responsible for lepton number conservation. In such a theory, the Majorana mass term 
is generated as the result of spontaneous breaking of lepton number. This theory can be the goal 
of future investigations; meanwhile, in the model under consideration, the asymmetry is believed 
to be generated at temperatures bigger than the electroweak symmetry scale, but lower than the 
scale where the Majorana mass is created. Thus the Lagrangian (|I]) is an "intermediate-energy" 
effective Lagrangian. At the high energies considered here, the vacuum expectation value for the 
Higgs condensate is very small, so that the masses (here we mean vacuum masses and neglect 
temperature contributions) of the charged leptons and Higgs particles are negligibly small 
or zero. 



2 Indirect contribution to leptogenesis. 

The Lagrangian (|Tp was considered several times 10] - 0], where the so-called "direct" contri- 
bution to leptogenesis was computed. Here we are interested in the "indirect" contribution to 
leptogenesis, calculated in references 0-||lO| and reviewed recently in Jll . 
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Figure 1: The self-energy diagram of the heavy Majorana neutrino 

The indirect contribution is described by the self-energy diagram in Fig.l. This diagram was 
investigated in the context of two different approaches. In the first one the self-energy diagram 
is considered to be the intermediate state of some physical process [|12|, |13] (for example, lepton- 



Higgs scattering [Tj]). In this case the four- momentum of the Majorana neutrino is off the mass 
shell and is determined by the four-momenta of the incoming (outgoing) lepton and Higgs boson. 
We wish to emphasize is that there are no restrictions on the value of neutrino four-momentum. 
In this approach the lepton asymmetry appears in the scattering process ~^ ^0^5 through 



interference of the tree-level diagram with the one-loop diagram |TT[ . 

The second approach, presented earlier in 0], uses the notions of transition Nj — > Ni, 
originating from the one-loop diagram mentioned already. If i = j then the diagram is analogous 
to any self-energy diagram and describes the correction to the fermion (Majorana neutrino in our 
case) mass. The precise role and calculation of this sort of diagrams is described in textbooks 



(see, for example, corrections to the electron mass in [|T^, plq] . 

A new situation appears for i ^ j. Now we have a free stable particle with mass Mj as 
incoming particle and it is impossible to transform it into a stable particle with mass Mj, because 
the process does not conserve energy and momentum. This transition can happen for unstable 
particles provided their mass difference is comparable to their widths. In addition the Majorana 
states are in a thermal bath interacting with the other fields so that energy and momentum is 
continuously exchanged with the background fields. If the transition probabilities are smooth 
functions of energies and widths, then it is a good approximation to calculate the transition 
amplitudes and probabilities per unit time and volume for unstable Majorana neutrinos and 
then introduce the results into the Boltzmann equations. In the latter step we also introduce 
the particle densities for initial and final states, as dictated by the thermodynamics of the early 
universe. A prototype of the thermodynamic calculation is described in ref. [Q, [|l7l and recent 



articles [p!8| , p!9[| . The interaction with the background fields is an integral part in the generation 
of the lepton asymmetry. In this article we calculate how the transition probabilities and the 
asymmetry are generated for unstable particles. 

Before we address this topic, we review how the mixing occurs in a few known cases. The mix- 
ing of the with is described by the box-diagram. In this case the i^'^-mesons are treated 
as composite particles made up of quarks. The mixing diagrams are computed at the quark 
level, which are the fields of the basic Lagrangian (Standard Model). In the diagrams appear 
vertices of the form usW, because the quarks mix in the Lagrangian through the Cabibbo- 
Kobayashi-Maskawa matrix. Thus the mixings in the above case originate from the mixings at 
vertices between quarks of different charges and manifest themselves as the mixing of mesonic 
states through the interaction of their constituent quarks. In the above case energy and momen- 
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turn is conserved. For the mesonic states the mass difference is comparable to the widths. For 
K-mesons 

1 1 

|AM^| = -\ATk\ Tk 




= 0.734 ±0.035. 



The neutrinos have charged current couphngs to the leptons where the lepton number is 
conserved to a large degree of accuracy. The Majorana couplings, on the other hand, may have 
lepton number violating charged couplings as seen in the Lagrangian of eqs. (1) or (3). The mass 
matrix of Majorana neutrinos in eq. (1) can be made real and at the classical (tree) level the 
generations do not mix. However, it was noted in ref. that the one-loop effects can mix the 
Majorana states, i.e. the states defined at tree level are not the physical ones. To state the result 
in another way, we find that Majorana neutrinos of different generation mix, while the electron 
and the muon do not mix. 

In this paper we show that the instability of the Majorana neutrinos is of principal importance 
for the mixing phenomenon. The neutrinos have finite widths and can be treated, through the 
uncertainty principle, as "fuzzy states" on their mass shell. This allows us to define the mixing 
of physical states in the case when the difference of the tree-level neutrino masses is less or about 
equal to their widths: 

\Mi - Mj\< max{Ti, Tj}. (5) 

Condition (5) can evidently be applied to any other particles. For electron and muon, however, 
rrin — irie ^ F^ and consequently the mixing is very small. 



and for the B-mesons |20 



3 S-matrix for the mixing of unstable particles. 

The general S-matrix theory necessarily uses the notion of asymptotic states. This means that 
states of incoming, \in{t —oo)), and outgoing, \out{t +oo)), particles are defined as states of 
free (noninteracting) particles at infinite times. The second principal element of the 5— matrix 
theory presupposes that interactions do not exist at all "infinite" times, but are "turned on 
adiabatically" from t —>■ — oo to t = and also "turn off" adiabatically from t = to t +oo. 
One, of course, understands that "infinite" here means macroscopically large in comparison to 
the time of interaction. This general argument can naturally be applied to the case of unstable 
particles. In this case the time of "turning on" and "off" the interaction is easily estimated to be 
of the order of the inverse particle width. We try to express this in a phenomenological definition 
(the idea is similar to that of Breit and Wigner when considering the cross section at resonance) 
of field operators for unstable particles. For neutrinos we introduce 



fc,A 



e'^'l*! (6) 



The indices i = 1,2 denote the generation of the neutrino and for the sake of simplicity we consider 
only two generations, a^^^ and a| ^ are the fermionic creation and annihilation operators. Formula 
(ID reflects the fact that neutrinos "disappear" (decay) for infinite time {Ni — > at t — >• ±oo). 
We assume that the widths are generated by the Lagrangian in eq. (1) and calculated by the 
decay diagrams or from the absorptive part of the diagonal terms of the self-energy. Eq. (6) is the 
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only step which is not proved from field theory, since field theory does not consider asymptotic 
states. We think, however, that the physical meaning of eq. (6) is rather clear and will help us 
understand the mixing of Majorana neutrinos. 

In this paper we will follow the idea of an earlier paper and will introduce an "effective 
Hamiltonian" with non-diagonal terms in the mass matrix. Initially they are calculated as 
elements of the S-matrix to second order in the Yukawa-couplings by using perturbation theory. 
Let i and j be neutrinos of definite fiavors, then 

Sij = J J d^x'{i\Lint{x)Lint{x')\j), (7) 
with the neutrino fields defined by eq. (H). 



4 Loop calculation with unstable Majorana neutrinos. 

As mentioned in the previous section, we shall use perturbation theory. To second order in the 
Yukawa couplings the S"— matrix element, corresponding to the one- loop diagram, is given by 
(1^. /^From the product of the interaction Lagrangians 

{i\ [hamlLa{x)(f){x)PRNm{x) + h*^^Nm{x)PL(j)^ {x)lLa{x)^ X 

X {h^Ji:{x')(j){x')PnN^{x') + hlj7:{x')PL<P+{x')lLa{x')} \j) 

two terms contribute to (|^): 

2hUhaj ■Wiix)PL(p+ix)lLaix) ■TZ{x')(Pix')PRN,{x') 

— — (8) 

+2h^^hl^ ■ PRN,{x)ct){x)lUx) ■ N,{x')PL(t>+{x')lU^'). 

The new aspect of our calculation is the form of the field operators for unstable neutrinos given in 
eq. (6). We calculate in detail the first term, using Wick's expansion for the product of operators 



-iq{x—x') 



(9j 



[ ik(x-x') ( 1 1 \ -i^p'2+M|t'-r, \t'\+ip'x'p j 

J (27r)4 [k^-ml + ie - + J 

where A is a parameter for the Pauli-Villars regularization.Q After integration over d^x, d^k, 
d^x' one obtains: 

c - \ [ ^+ [ rl+'oh* h P f /j + ^a f dk^ \ 

with a similar expression for the terms with the Pauli-Villars regularization parameter. An 
important feature of this formula is the product of the two integrals 



^Another regularization is equally possible. 
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In the case of vanishing widths Fj = F,- = each of the integrals is equal to a 5— function 



and their product S{^Jp'^ + — — q^)5{^'p'^ + M? — kP — g°) is zero except when Mj = Mj 
(with f) = p'). This result is a mathematical demonstration of physical arguments, discussed in 
section 2. On the other hand, for non-zero widths the integrals Ji and J2 lead to "bell-shaped" 
Lorentzian functions: 

Ji • = {2nr ^^it . (12) 

Ff + {^p + Mf - F - g0)2 F| + {^p + M] - A;" - q^Y 

This product is non-zero and large when condition (|^) is satisfied; its non-zero value can be 
graphically understood as the "overlap" of two "bells". So the limit |Mj — Mj\ S> Fj(j), which is 
often used in the calculation of asymmetry, leads to very small overlap functions. Another remark 
concerns the three-dimensional 5— function with zero argument (27r)^5(p— p') = (27r)'^5(0), which 
appears in eq. (10). We will keep it throughout the section and it corresponds to the volume 
element where the interaction takes place. 

When we substitute (0) into (p!0D, and in addition integrate over we obtain the first order 
expression for Fj(j): 




■X 



P + Mf- + M]f + F 
1 



-g0)2-(p-g)2_^2 (^p*2 + M2 - g0)2 - {p- q)^-h?) (13) 
Vi 1 



^ {^p + M] - y'p2 + jvf2)2 + r2 

1 1 



X 



IP + M] - g0)2 - {p- g)2 _ ^2 (^^2 + JVf2 - g0)2 _ ^2 _ A2^ 

This expression, as is easily seen, reproduces the usual expression for the loop integral after one 
introduces the zero- component of four momentum p^: 



(27r)4 q^-ml + ie " 




1 1 



p2 + M2)2 + F2 (p° - q^f -ip-q) 
1 1 



2 — m'i 



(14) 



Alf 



(/ - + Mff + F2 (/ - - {P- q)^-rr 
Integrals of this type are standard and the result for the limit 0, — is well known: 

d^q /i + nia f I 1 \ 



Int = lim 



mc,m^^O J (27r)^ g2 _ 777,2 -|- ^5 ^ (pM — gM)2 _ ^772 _(_ ^^^i _ g,/^^ _ y\2 _|_ 

9dis - 7;9abs , 9dis = -77-^ o ^ + 7 ' S-afes 



2"""V' " 167r2 \ 2 p2 '4;' istt 
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It is recognized that the dispersive part of this integral (gdis) can be reabsorbed in the definition of 
couphng constants, while only the absorptive part (gabs) survives and has physical consequences; 
for more details see |]15|, |l^] and references therein. 
Now the matrix element Sij is given by 



Sij = hUhaj27i6{p-p){27rf{-i)x 



-i/2- gabs 



TT 



IP 



p + Mf- + M2)2 + r 

■ gabs 



j 



We transform this expression with the help of the Dirac equation 



pa 



uI^PrM,. 



7Vp' + M]-^P Pru^^, = M,PluI„ 



vL,Pl 



pa' 



7%/p2 + M2 - 7p 



-v'pa'PRMj, 



Finally we arrive at 



sip = hi^K,27,5{p~m^)\-^rA { 



^yp + M^-ip 



1 , 



PrV;„ = -MiPLU 



pa' 



pa 



-i/2 ■ gabs 



TC 



+ M, 



-i/2 ■ gabs 



' (. Ip + - Jp + Mff + Tj ' 



(15) 



(16) 



;i7) 



This is the final expression for the S-matrix element originating from the first term in eq. 
The second term from eq. (|^) leads us to the similar expression with f — spinors 



-i/2- gabs 



^ {^p + - Jp + M]Y + V] 



-i/2- gabs 



^ up + M2 - Jp + Mjf + 



pa 



As mentioned already, (5(p — p)(27r)^ represents the volume element V . Similarly the time of 
interaction also appears as multiplicative factor. 



271- 



(r. + r,)/7r 



(r, + r,)2 + {^p + - Jp + M]Y 



(19) 
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The new feature of this calculation is the presence of the masses and the decay widths in the 
S-matrix. The new terms define the time intervals when the interaction takes place. To obtain 
the transition amplitudes per unit volume and unit time we must divide by the factor V ■ Jij 
where V is the volume element and Jij the time interval of the interaction. In other words, we 
shall work with the "T-Matrix" defined through the equation 



where the energy conserving (5-function will be substituted by the expression 

2(r. + r,) 



(20) 



(21) 



One easily arrives at 

oil) n{n) 



(-1) 



2 16 



— (h* h 

67r \"'ai"' 



Mi ^ „ Pr + M, „ Pr 



M. 



u 



PR + Mi 



Pt 



pa' 



^pcr 



(22) 



Notice the S-matrix element (17), (18) vanish in the case |Mj — Mj\ 3> since the Lorentzian 
representation of the Delta-function (21) tends to zero. It is easy to see that the same occurs 
for the transition probability 



V 



Jij 



IT- 



(2vr)^-5(p-p') 



(r, + r,)/vr 



iEj-E,)^ + iTj + T, 



Let us summarize our results. We have shown that neutrinos of different generations can mix 
through one-loop diagrams even if they have different masses. The physical reason of this 
possibility is related to the final widths of neutrinos : their mass shells are "fuzzy" and energy 
is conserved with the accuracy comparable to their width. The closer the masses of neutrinos, 
the larger is the transition Ni Nj probability; and for |Mj — Mj\ ^ ^nj) mixing has no 
physical meaning. 

As we already know, neutrino mixing provides an indirect contribution to leptogenesis. As 
it was previously done in f^, 0, we are working in terms of an effective Hamiltonian, or, what 
is the same, in terms of an effective mass matrix. From the transition matrix elements and the 
mass term in eq. (1) we can now proceed to calculate an effective mass matrix and the physical 
states. 

We consider the expression for Tij as the first order contribution for mass matrix elements. 
One should notice that the contributions are different for R and L parts of our 4-spinor Ni in 
eq. (2), i.e. they are different for A^'^ and {Nm)'^. Evidently, an overall factor of (— i) and the 
spinors originate from the definition of the S-matrix and do not appear on the mass matrix. A 
factor of 1/2 is also omitted in order to be consistent with the Lagrangian in eq. (1), which has 
1/2 in front of the mass term. 

So the corrections to the Mi [NjnY Njn mass term (see Lagrangian in the form (3)) are 



Hij 



2 8 



TC 
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and the corrections to the Mi Nm {Nf^iY mass term are 



2 svr 1, r, + ^ r, + r, J 

When Fj = Tj, this resuh coincides with those obtained in ref. 0. Our account of finite widths 
of Majorana neutrinos shghtly corrects the result which is not very essentiaL What is important, 
is that the definition in eq. @ enables one to give physical meaning to the mixing of particles 
with different masses. 



5 Effective Contribution to the Mass Term and Physical 
Neutrino States. 

We shall work in a vector space with four basis vectors (iV^]^, Njn, -^ra) which are the 
states occurring in the interaction term of eq. (1). The transitions among these states introduce 
an effective matrix 



m 






Ml + ifii 





H12 


Ml + Hu 


















M2 + H, 


Hu 





Ma + H22 






(23) 



with 



h*.h 



+ h h ^ 

\ I I'm' I' n 



-HI 



and 



1 



9ab 



9ab 



IGtt 



(24) 
(25) 



The 2x2 matrices in the upper left and lower right corners of eq. (23) are lepton-number 
violating but flavor conserving, while the 2x2 matrices along the off-diagonal are lepton- and 
flavor-number violating. All terms should be present and we found no way to reduce it to a 
2x2 matrix. It remains to diagonalize the matrix and flnd the wavefunctions. It is instructive 
to present a perturbative solution of the problem and then discuss the exact solution. 
We split the mass matrix into the dominant term and a perturbation 



m 



H. + XV 



with 



Hn 



/ Mi + Hu 
Mi + Hn 

M2 + H22 

V M2 + H22 



and 



/ 



XV 







-XH 



12 






XH12 





-XH 





12 



XH^2 \ 



/ 



(26) 



(27) 
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We introduce a small parameter A in order to keep track of the perturbative corrections and 
at the very end set A = 1. Applying perturbation theory we find the eigenfunctions to O(A^) 



I 1 \ 

1 

X 
\ Y ) 



with the eigenvalues Ai 



Ml 



2 ' 



A, 



and 



Mo 



,•£22 



( X'\ 
Y' 
1 

V 1 / 



and Rii 



(28) 



jFij. They have the 



time dependence f/j e * 2")*^ with % = \ and 2 which is consistent with the time dependence 
introduced in eq. (6). The physical states are 



^1 



\Nnx)^Wm)^X\N%,)^Y\Nni)\ 



^2 = -^[|iV^2) + |iV^2)+^'|iVm) + >^'|A^m)] 



(29) 
(30) 



with X = , Y = and similar formulas for X' and Y'. 



M'(-M'^ 



The results are the same as in ref. [5, 6] with an additional dependence on the widths. The 
interesting result is that the definition in eq. (6) makes it possible to give a physical meaning to 
the mixing of particles with different masses. 

We can finally calculate the asymmetry for the decay of each state. The width for the decay 
of ipi into leptons is 

r(^, ^£+...)oc^|/i„i + /i„,F|2 (31) 



and into antileptons 

r(^i^^+ 

The lepton asymmetry, defined as 
and it is straightforward to calculate it 



oc 



(32) 



(33) 



5i = 7^ 



Ml M2 



Im (/i*i h 



a2 



Stt Ml -Ml \Ki\'' + \K2X\'' + ReK^K2{Y + X* 



(34) 



In the above calculation we considered the case IM2 — Mi| ^ H12. In the case that the two 
Majorana neutrinos are nearly degenerate we must diagonalize exactly the matrix in eq. (23) 
and then we recover the resonance phenomenon introduced in eq. (17) of ref. 0. 

Our study so far considered the neutrinos as free particles. In reality, however, they are in 
a backgound of fields interacting many times with the other particles. These interactions at a 
finite temperature can modify their masses and widths []T8|, |19|. In the present analysis we have 
found that within the range of validity of eqs. (18) and (19) the T-matrix in eq. (22) is a smooth 
function of masses and the widths, and for this reason we can introduce the matrix elements in 
the Boltzmann equations in order to study the development of the asymmetry. This approach 



was followed in ref. \\A . 
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Summary 



We have shown that finite widths of Majorana neutrinos play a principal role in producing 
neutrino mixing via self-energy diagrams. In our approach the widths of neutrinos are treated 
as "fuzzy" states on their mass shell and the mixing is understood as an "overlap" of two fuzzy 
mass shell states. For the mathematical reahzation of these ideas we changed phenomenologically 
the definition of asymptotic neutrino states, including their widths as in eq. (6). Concerning 
the results, we showed that these changes do not infiuence significantly the expressions for the 
effective mass matrix and for the asymmetry generated in the case of small mass difference. For 
large mass differences the asymmetry is suppressed. 
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